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Abstract. We analyze a system of interacting islands of XY spins on a triangular lattice. This model has 
been introduced a few years ago by Eley et al. to account for the phenomenology in experiments on tun¬ 
able arrays of proximity coupled long superconductor-normal metal-superconductor junctions. The main 
features of the model are the separation of a local and a global interaction energy scale and the mesoscopic 
character of the spin islands. Upon lowering the temperature the model undergoes two crossovers corre¬ 
sponding to an increasing phase coherence on a single island and to the onset of global coherence across 
the array; the latter is a thermodynamical phase transition in the Ising universality class. The dependence 
of the second transition on the island edge-to-edge spacing is related to the proximity-effect of the coupling 
constant. 


1 Introduction 

Recently Eley et al. [I] have introduced a model of cou¬ 
pled islands of XY spins. Their goal was explaining the 
results of measurements of resistance in arrays of long 
superconducting-normal-superconducting junctions. The 
experimental devices are based on planar arrays of iden¬ 
tical islands made of superconducting (Nb) grains, dis¬ 
posed in a triangular matrix over a metal (Au) film. The 
authors studied the dependence of the system resistance 
R{T ; h, €) on temperature T, on island (vertical) thick¬ 
ness h and on inter-island spacing i. They found i) the 
resistance dropping to zero, by lowering T, in two steps, 
and they determined two transition temperatures Ti and 
T 2 with Ti > T 2 ; ii) an interesting dependence of both 
Ti and T 2 on the island spacing, possibly (both) extrap¬ 
olating to T = 0 at large iii) a strong dependence of 
the behavior of the system on the island thickness. In a 
following paper [2], they discussed a more detailed com¬ 
parison between the experimental data and the predictions 
about the dependence of T 2 on I given by the conventional 
theory of Lobb, Abraham, and Tinkham (LAT) [3]. They 
argued that for large inter-island spacing the supercon¬ 
ducting transition is more likely to be driven by diffusive 
effects PIIS] in the normal metal substrate, and that it 
does not depend on the details of the superconducting is¬ 
lands, with the puzzling dependence on island height as a 
notable exception. 

The superconducting transition in proximity-coupled 
macroscopic grains embedded in normal metal films has 
been the object of intensive work in the past years, [51171 


IH]- Tunable realizations of 2D superconductivity were also 
object of previous experiments [9]. The classical model 
presented in Ref. [I] to account for a novel phenomenol¬ 
ogy is at difference with previous theoretical and experi¬ 
mental work, as it takes into account the intrinsic fluctua¬ 
tions of the superconducting state inside the single meso¬ 
scopic islands (see also refs. mm for recent theoretical 
and experimental work on mesoscopic Sn islands laid on 
graphene). It is clear that, because of many reasons we 
will discuss in the following, this model does not try to re¬ 
produce faithfully the experimental situation (for example 
the use of an anisotropic coupling is not connected to the 
physical form of the Josephson interaction but is a tool 
needed to obtain a phase transition). The idea of [T], and 
our point of view here, is to analyze a very simple model 
that offers a behavior quite similar to the one detected 
in the experiments, and to try to learn from this behav¬ 
ior. Here we will present a detailed analysis of the model, 
that corroborates and supplements the hints coming from 
the first analysis of [I]. It is also worth mentioning that 
tunable two-dimensional superconductors are also of inter¬ 
est in a revived search for a non-conventional 2D, T = 0 
metallic phase. [niiKiKisEaiiziiiaiis] 


2 The model 


The Hamiltonian of the model is based on 0(2) vectors 
living on the individual grains (labeled with i,j). Groups 
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of grains form islands (labeled by p): 


h = -jY1 E 

p {ij)&p 


- Mp- yMp, , 

(1) 

{p,p') 


III 

M 

Cq 

(2) 

iGp 



where by a dot we denote the scalar product in the internal 
space and where J' is a 2 x 2 matrix of couplings. Each 
island is a D-dimensional hyper-cubic array of grains of 
linear size I (and volume V/ = , with either D = 1 or 

D = 2 in our computations). Islands are arranged on a 
(two dimensional) planar regular lattice of linear size L. 
Islands are mesoscopic: their linear size / is not larger than 
a few grains. Because of that they may have large global 
phase fluctuations. The size of the underlying planar array 
is macroscopic, L ^ I. The case of one-dimensional island 
is an exercise useful to understand better the role of island 
dimensionality, and does not try to be a description of the 
experimental situation. On the contrary the case of two 
dimensional islands is probably closer to the experimental 
situation, where islands have many layers, but only one or 
few conjure to build the inter-island interaction. 

The first term of the Hamiltonian is a sum of nearest- 
neighbor interactions between grains contained in the same 
island. The second term couples neighboring islands in the 
array. Each spin in a given island interacts directly with its 
neighboring spins in the same island and with the average 
spin field of surrounding islands. In the model proposed 
in Ref. [T] the inter-island coupling matrix (in the internal 
0(2) space) J' is anisotropic: 



This particular choice polarizes the islands in one spe¬ 
cific direction in internal vector space, changing the na¬ 
ture of the inter-island phase transition. This is a techni¬ 
cally useful choice (since it carries a phase transition in 
the game), but it does not aim at reproducing the de¬ 
tails of the physical Josephson interaction. Finally, notice 
that in the isotropic case and if the energy scales J and 
J' are far apart, i.e. J' <C J, so that at low tempera¬ 
tures all (mesoscopic) islands are magnetized, we recover 
a Kosterlitz-Thouless [20] phase transition. 

The island-island couplings depend on the temperature 
and on the inter-island edge-to-edge spacing, according to 
the theory of diffusion of electron pairs in SC-Normal-SC 
junction. As in the work of [T] we take a “quasi-proxfmity- 
effecV' [213] form for both couplings; in a proximity in¬ 
teraction J' would depend on the inverse square of island 
spacing when the latter is small, but following [2 and 
for the same sake of simplicity we omit this part of the 
interaction, that is not expected to change the nature of 
the phase transitions here. We assume the proximity-effect 
form also for the grain-grain coupling J and we take the 


grain-grain distance as the length (lattice) unit (and de¬ 
note the inter-island spacing as €}. 


J = Jo exp ! 

(4) 

J' = Jo exp . 

(5) 


The choice of an interaction of a proximity-like form im¬ 
plies that physically grains of the islands are also im¬ 
mersed into a metallic matrix. The authors of Ref. [2 in¬ 
troduced the model defined in o to explain the presence 
of two transitions (intra-island, Ti, and inter-island co¬ 
herence, T 2 ) and the depression of Ti for increasing island 
spacing 1. Such a dependence of Ti on I has been observed 
and reported for the first time in [2 (for example in previ¬ 
ous experiments on lead disks on a thin substrate [3] where 
islands were not mesoscopic, the effect was not observed). 
The energy scales Jq and Jq must be well-separated: we 
adjust them in order to clearly split the high-T and the 
low-T transition. We fix Jq = 1 and vary Jq; in order to 
easily compare data for different island sizes, we also take 
Jq = Jo/Vr and adjust the parameter Jq. 

In Ref. [2 the authors also give some predictions by 
analyzing a D = 1 islands model, where it turns out that: 

— T 2 < Ti provided J :$> J' and islands are small; 

— Ti —0 when J' —>• 0. 

The second statement is rather counter-intuitive, as it im¬ 
plies that islands are not superconducting when they are 
isolated. This implies that an array of superconducting is¬ 
lands can be superconducting even if the inter-island spac¬ 
ing is larger than the superconductor coherence length, 
but an isolated island of superconducting grains, where 
grains are packed closer than islands are in the array, 
looses phase coherence. In this respect, the one-dimensional 
and mesoscopic character of the islands plays a role, since 
for macroscopic chains we must expect Ti ~ 0, and, as 
we will see in the following, for large / the intra-island 
coherence is driven by inter-island ordering (also see the 
discussion in the Conclusions section). 

Another striking aspect of the phenomenology of the 
system[2|2] is the dependence of its behavior from the 
height of columnar grains. Realistic islands extend in more 
than one dimension. We have analyzed by numerical sim¬ 
ulations the behavior of one and two-dimensional islands. 
The dependence on thickness may suggest that it would 
be interesting to go to D = 3, too (in case of a very large 
value of / this should turn the Ti transition to a true 
second-order one). If energy scales are adequately sepa¬ 
rated (i.e. J » J'), this should not change the properties 
of the T 2 (KT) transition, when phases of grains in the 
same island are mutually locked. Mesoscopic islands can 
then have a crossover at Ti from a disordered to an or¬ 
dered phase; for D > 2 and large I this crossover is related 
to a true thermodynamic transition. 

3 Numerical simulations 

We have studied the model defined in ((2). By following a 
very slow annealing protocol, with constant ratios between 
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adjacent temperature (a logarithmic annealing scale), we 
have cooled down the system in order to get a signal for the 
two transitions. At each temperature we collected mea¬ 
surements during the evolution of the Monte Carlo dy¬ 
namics. Our Monte Carlo step consists of sweeps of the 
whole lattice by single-spin moves Metropolis dynamics, 
followed by Uo sweeps by over-relaxation [22]. The choice 
of Um = 10 and Uo = 12 have shown to be appropriate 
for most island and array sizes considered (and an overkill 
for the smaller sizes), and allowed an estimate of inte¬ 
grated auto-correlation times not larger than ten Monte 
Carlo steps at most temperatures. Although averages al¬ 
ways stabilize quickly after any temperature change, we 
drop the first half of the collected measurements at all T 
values. The simulated annealing protocol, together with 
over-relaxation, is appropriate to the needs of this prob¬ 
lem. All observables of interest converge very fast to a 
plateau at all temperatures. Although averages always sta¬ 
bilize quickly after any temperature change, we drop the 
first half of the collected measurements at all T values. 

Since the devices in the experimental setup m are 
triangular arrays, we consider a triangular lattice. A sim¬ 


ple implementation choice in simulation is to consider a 
triangular array with regular hexagonal shape with heli¬ 
cal boundary condition (in this way we preserve the sym¬ 
metries of the triangular array and avoid involved bulk 
properties extrapolations); each side of the hexagon has a 
width of L islands, and the number of islands is V 5 = 
3L{L — 1) -I- 1. We simulated systems with L = 8,16 
and 32: for I? = 1 systems we have islands of sizes I = 
16,36,64,100 and 144, while for Zl = 2 we have I = 
6 , 8 ,10 and 12. The inter-island edge-to-edge spacing i 
has been varied in the set {2,4,8,12} for D = 1 and 
{2,4,8,12,16,24} for £> = 2 (the SNS arrays in Ref. [T] 
had edge-to-edge spacings up to approximately 10 times 
the grain size in their experiments, and £ up to 20 in 
Ref. 0)- We have considered both free (FBC) and pe¬ 
riodic (PBC) boundary conditions on the single islands. 
Although we found no qualitative differences, FBC is a 
more realistic choice when dealing with mesoscopic ob¬ 
jects, for which we expect finite-size effects to play a role. 
We measured the following quantities. 

— The single island magnetization magnitude (averaged 
over islands): 



Mi 



Vr 


i€p 


( 6 ) 


This should be, in the limit of infinitely extended is¬ 
lands, a good order parameter for island internal order¬ 
ing (in any direction in internal spin space, and glob¬ 
ally over the array: it has a non-zero value whenever 
any island starts to order internally and it is maximum 
when all islands are locally ordered, independently of 
the relative orientation between different islands). 

— The total magnetization: 


M = 


VsVi 




p i€p 


(7) 



— A renormalized magnetization: 


— 



( 8 ) 


which is a unit vector on the single island, and its 
average over the array. 

— Af/j, which characterizes the globally-ordered phase, 
even if islands are not yet internally fully ordered: 

Mr = E 


We also consider the fluctuations of the magnetizations 


Fig. 1 . D = 1 results for the largest simulated size L = 32. We 
plot the connected susceptibilities xi s-nd xr for t^'o different 
values of I. jo = 0.0072. On each panel we also plot the curves 
for four values of the island spacing £. Notice that the peaks 
of xi (which mark the Ti transition) depend strongly on £. 


X = ViVs [{M^} - (|M|) 
XI = Vi [{Mf) - {Mif 
XR = Vs\{Mli)-{MRf 


( 10 ) 

( 11 ) 

( 12 ) 
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where (M|) and {Mj) are averaged over all islands, x is 
the total susceptibility of the system. At very low temper¬ 
atures, when Ml ~ 1 , we have XR ~ xl^i- At T 2 , that 
we define as the location of the peak of the inter-island 
susceptibility XRi XR y have very similar sharp peaks 
(both in shape and location). We take the location of the 
(very smooth) maximum of xi ^s a rough estimate of the 
temperature Ti at which islands order internally (in this 
way we give an operative definition of Ti in our model: 
since islands are of finite extent the Ti defined in this way 
is indeed a crossover temperature). 


4 Results and discussion 

In Fig. [T] we report the results for arrays of one-dimensional 
chains. The Ti temperature value goes to zero very fast as 
the island size grows, as expected for linear spin chains. 
Upon lowering T, coherence between island builds up and 
also drives the internal ordering; the two transitions can 
be resolved only for very small island sizes and by lowering 
considerably the value of coupling constant The effect 
is also strongly dependent on island size. 

The situation is far clearer for two-dimensional islands 
(see Fig. El), where we still have a finite temperature ther¬ 
modynamic transition for isolated islands in the limit of 
large sizes. For mesoscopic islands, the crossover between 
unordered and ordered island depends more weakly on is¬ 
land size than in the linear chains case. Our numerical sim¬ 
ulations show that the temperature Ti does not depend on 
the island spacing, or the dependence is very weak. This 
effect has been also reported in experimental results on 
non-mesoscopic island samples [9]. Also the dependence 
of Ti on island size is very weak. 

We try a more quantitative approach studying the de¬ 
pression of T 2 by increasing the inter-island spacing. We 
take as an estimate for T 2 the midpoint of the temperature 
interval bracketing the peak at its half-height. The depen¬ 
dence of T 2 on i and I for the largest simulated array size 
L = 32 is shown in Fig. [31 

Following Ref. [T], we notice that T 2 {£) compares well 
to a proximity-effect prediction 

T 2 = Z\exp(-C£v^) , (13) 

corresponding to the solid straight line in Fig. |3l suggest¬ 
ing a diverging £{T 2 = 0); we report in Table [Hour best 
fit estimates of the parameters A and C. 

We have run more accurate numerical simulations in 
the temperature region close to the T 2 transitions, with a 
four times smaller cooling rate and ten times more mea¬ 
surements. We measured the Binder cumulant 



The value of G 4 at the T 2 transition point is universal [22] ; 
we report data for / = 6 , £ = 4 and various array sizes L 
in Fig. [H Note that the dips in the curves of G 4 for the 





Fig. 2 . D = 2 results for the largest simulated size L = 32. 
We plot the connected susceptibilities xi (filled symbols) and 
Xr (empty symbols) for three different pairs of the parame¬ 
ters I and Jq. On each panel we plot the curves for six values 
of the island spacing 1. The peaks of xr (which mark the Ti 
transition) depend only weakly on i. In all panels J'q = 0.072 


largest system sizes in Fig. [Tj are due to the breakdown of 
the 0 ( 2 ) internal symmetry introduced by the anisotropic 
form of J' in Eq. ©■ This is the behavior one would ex¬ 
pect, since at high temperature the G 4 value depends on 
the symmetry of the system. The fluctuations of the mag¬ 
netization in the infinite volume limit are Gaussian dis- 
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numerical evidence for a second order phase transition in 
the two-dimensional Ising universality class. The asymp¬ 
totic value of the crossing points of the Binder cumulants 
curves (see inset of Fig. 0 , T 2 c(T, 2L) which asymptoti¬ 
cally tends to T 2 , is clearly different from zero. We finally 
remark that the value of the Binder cumulant below the 
critical temperature is almost unity as expected asymp¬ 
totically (as the size of the system goes to infinity). 


Fig. 3. A scaling plot of the T 2 transitioir temperature for 
arrays of D=2 island based on the behavior Aexp{—C£T 2 ^^) 
as suggested by Eq. [5] (L = 32 data). 


tributed and the Binder parameter for the two-component 
magnetization of a XY system should approach the value 
G4{T —>• 00) = 0.5, whereas for Ising spins the correspond¬ 
ing high-temperature value is G 4 (T —>■ 00) = 0. When long 
range order in the system builds up at low temperatures, 
the value of the Binder parameter must approach unity: 
G 4 {T —>■ 0) = 1. The data in Fig. |3]show that G 4 is not 
monotonically increasing when the temperature decreases: 
it starts at a value around 0.5 but, as soon as the inter¬ 
island term becomes more important with respect to the 
intra-island interaction in the Hamiltonian, the effects of 
the Ising symmetry sets in and in proximity of the critical 
region, just above T 2 , the value of G 4 drops to low values, 
as expected for an Ising system. The minimum of the dip 
decreases as the system size L increases. 

Moreover, the critical value of the Binder cumulant 
(which is universal) for the two-dimensional Ising model is 
known to great accuracy 



0.32 


0.36 


0.4 


0.44 


Fig. 4. Binder cumulant G 4 versus T for D = 2-islands with 
1 = 6 and £ = 4 and for the five simulated sizes (L). In the inset 
we show in detail the region near T ~ 0.335, where different 
Binder curves cross. 


5 Conclusions 


It has been very difficult to resolve the crossover at Ti 
(internal island ordering) and the T 2 transition (inter¬ 
island ordering) in the case of one-dimensional islands. 
. The values G 4 = 0.9160386(24)In D = 1 and in the limit of large islands sizes we ex¬ 


compares extremely well with our value of the Binder pa¬ 
rameter at crossing, close to T 2 ~ 0.335 (for comparison, 
from the position and width at half-height of the peak of 
the susceptibility for the same simulated system we ob¬ 
tain T 2 ~ 0.340 ± 0.004) We report in Fig. [5] the details 
of the crossing of the Binder curves. This provides a clear 


I 

A 

C 

6 

3.74(10) 

0.935(6) 

8 

5.21(15) 

0.890(6) 

10 

7.56(28) 

0.889(7) 

12 

8.15(26) 

0.840(6) 


Table 1. Best fit estimates of Z\, G parameters in Eq.[5]from T 2 
data for various island size I (D = 2, L = 32, £ = 8,12,16, 24). 
Data points for the shortest distances £ = 2,4 have been ex¬ 
cluded from the fits. The chi-square per degree-of-freedom val¬ 
ues vary between 2.4 and 3.7 and the quality-of-fit parameters 
between 0.11 and 0.32. Uncertainties are gnuplot estimates cor¬ 
rected as in m- 


pect Ti —^ 0 and no thermodynamic transition at finite 
temperatures. In the range of simulated island sizes the 
measured crossover temperature at which the mesoscopic 
island order internally is as small as T 2 . The island inter¬ 
nal magnetization remains small and no clear maximum 
of the susceptibilities signals a crossover down to T 2 . At 
T 2 the inter-island interaction couples the fluctuations of 
the magnetizations of neighboring islands, making them 
coherent: at this point spins inside each islands starts to 
align to the average field of neighboring islands. The meso¬ 
scopic character of the islands is crucial: in the limit of 
large islands we expect the fluctuations of local magneti¬ 
zation to be too small (we did not try an experiment in 
that direction). Then, although J' <g; J, it is the T 2 tran¬ 
sition that drives both inter-island and internal ordering. 
As the inter-island spacing grows the depression of T 2 im¬ 
plies the depression of Ti, too. This is compatible with the 
counterintuitive requisite that Ti —>■ 0 as J' —>■ 0 discussed 
above. 

Since islands in experimental setups are not laid down 
on the substrate as unidimensional chains of columnar 
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0 0.12 r™ 0.24 0.36 


Fig. 5. Binder cumulant G/i{T2c{L,2L)), for D = 2-islands 
with 7 = 6 and 7 = 4, computed at the crossing point 
of the Binder curves of lattice sizes L and 2L, denoted as 
T2c{L,2L), as a function of We have marked, with a 

green horizontal line, the Gi value for the two dimensional 
Ising model: G 4 = 0.9160386. In the inset we show the behav¬ 
ior of the temperature-value of the crossing point T 2 c{L, 2L) of 
the Binder curves as a function of We have assumed 

the scaling of the two dimensional Ising model: T2c{L,2L) — 
T 2 ~ and G4{T2c{L, 2L)) - Gi{T 2 ) ~ We have 

used the exact value v — 1 and the conjectured one for the 
correction-to-scaling exponent A = 7/4. For more details, see 
Ref. [23]. 


grains, the ID system is not really connected to the exper¬ 
imental setup we analyze here. The 2D island system is, 
on the contrary, closer to the physical system we want to 
understand, and resolving the two transitions is easier in 
the case of the D = 2 system. The Ti transition for D = 2 
macroscopic islands and for J' J is expected to be of the 
Kosterlitz-Thouless type. Since we consider mesoscopic, 
finite islands, we observe, as expected, a “long-range” or¬ 
dering of the islands. We expect a spin ordering crossover 
at a temperature value Ti which does not go to zero with 
the island size as fast as in the D = 1 case. The model 
with planar islands is capable of describing the depression 
of the global superconductivity transition down to very 
low temperature at large lattice spacing. The transition 
temperature we measure by locating the peaks of response 
functions remains finite for moderate-to-large inter-island 
edge-to-edge spacing. T 2 approaches a zero value only for 
very large inter-island spacing. We cannot detect any de¬ 
pendence of the Ti transition temperature on At small i 
values, the value of Ti is influenced by nearby islands only 
because of global coherence driving internal ordering. We 
did not try a full finite-size scaling analysis of the T 2 tran¬ 
sition, which for well-separated coupling scales jo and jg 
is in the Ising universality class. On the simulated length 
scales there are no relevant effects of the mesoscopic char¬ 
acter of the islands on the behavior of which behaves 
at large inter-island spacings as expected by the choice of 
the dependence of the coupling constants on i and T. 


Appropriate variations of the basic model we have dis¬ 
cussed here could lead to interesting developments in the 
study of the superconducting transition in arrays of SNS 
junctions. We think it is an interesting starting point to 
understand many striking experimental evidences, as for 
instance the dependence of the transition temperatures on 
island thickness, or the strong depression of the Ti and T 2 
transitions. 

We thank Kay Kirkpatrick for introducing us to the model 
studied in this work and her and Jack Weinstein for interesting 
discussions. This work was partially supported by European 
Union through Grant No. PIRSES-GA-2011-295302, and ERG 
Grant No. 247328, by the Ministerio de Giencia y Tecnologia 
(Spain) through Grant No. FIS2013-42840-P, and by the Junta 
de Extremadura (Spain) through Grant No. GRU10158 (par¬ 
tially founded by FEDER). 
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